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ABSTRACT 


Three regression techniques for the calculation of ballistic winds 
using satellite infrared spectrophotometer (SIRS) data are evaluated. 
Experiments testing the derived methods are carried out using two 
SIRS-B data set. In all three cases, the data sets are divided 
into a large dependent sample used to derive the regression equations 
and smaller independent sample for testing the equations! validity. 
One method, that of geostrophic calculations using regression 
derived "D" values of ballistic geopotential height is found to be 


superior to the other methods. 
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I. INTRODUCTION 


The resultant wind which deflects a reentry vehicle as it 
travels along its flight path is known as the ballistic wind. This 
is a weighted integral of the winds encountered at the various levels 
of the atmosphere through which the projectile travels. The general 


momo. ohne ballistic wind calculation is: 


Poo 


Wrine W (inp) W (In p) d(\n p) 





where W,, is the ballistic wind, p is pressure, Poo is surface pressure, 
WV is the actual wind at level p, and ee is the ballistic wind 
weighting factor for a Layer. ME is proportional to the effect of 

the winds in a particular layer and the time the vehicle remains in 
that layer. A similarly computed quantity, ballistic density, is 
representative of the total retarding force of the atmosphere on a 
reentry vehicle. 

Present methods of evaluation of ballistic winds and density 
involve the use of rawinsonde data and rocketsonde data as available 
for the levels in the integration of Eq. (1). Smith, et. al. [6} 
have recently reported success in obtaining atmospheric temperature- 
pressure profiles using satellite infrared spectrophotometer (SIRS) 
data. This success implies that ballistic winds and densities could 
be calculated in any region without the necessity of a rawinsonde 


release. 





Ruggles [5] has noted that the envelope Be ieriine functions 
used with the individual SIRS channel radiances in obtaining 
atmospheric soundings and the ballistic wind weighting factors are 
analogous. He further suggested that the ballistic winds and 
densities might be calculated directly from SIRS data without first 
calculating the temperature-pressure profile. Ruggles! development 
is outlined below. 

It is the object of this research to investigate the feasibility 


of Ruggles! suggestion with respect to ballistic winds. 





TIT, DEVELOPMENT OF EQUATIONS 


A. FHASIBILITY OF THE USE OF SIRS DATA 
If the geostrophic approximation is assumed, Eq. (1) can then 


be expressed as: 


Kin (3 —=— {k x Me 2(x x) 4) Inp) \W / (np) d((n °) (2) 


or, rearranging terms: 


== ae (2 4) inp) W, (inp) d Cin go) ©) 


F Inp. 


where g is the acceleration due to gravity, f is the Coriolis 
parameter, Re is the horizontal gradient of geopotential height 
on a constant pressure surface. Now, expressing the weighting 


function W as a linear combination of known arbitrary functions 


of inp, F, (1n p): 


W (in p) = z a. F. (in p) (4) 


where a. are constants. Similarly, Z(x, y, ln p) can be expressed 


as ¢ 


n 
Z(x, y, 1n p) = S ol He Ge y, in p) (5) 
i 


where r. are the constants and He are a set of arbitrary functions 


Ciinepe Substicuiion of lq. (4) and Eq. (5) into Faq. (3) produces: 








MF, ke £:r; [7:64 ne Fi (np) dnp) (6) 


Now consider the radiative transfer equation for the radiance 


observed by the ih channel of the DIRS, sensing at a wave number, V ‘ 


A(ing) (7) 





N.= No . [A,V, Ting) TO, Inv] {8 Ly, T(Inp)] 


3 2 inp 
In p) 


Here, N; is the radiance sensed, and oe is the radiance emitted 
by a contaminant such as the top surface of clouds. This “lo is a 
function of A, the amount of cloud cover, y » the wave number of the 
it channel, and T(1n p) the temperature structure above the clouds. 
The Planck radiance at wave number, » and temperature, T, is 
represented by By LU T(1n p)} - and{ (y, In p) is the fractional 
transmittance of the atmosphere at wave number y from pressure level 
p to the SIRS. 

Smith et. al. [6] have shown that cloud contaminant term can be 


evaluated by an iterative procedure. With that term evaluated, 


(7) becomes: 


-[ mi lV sTCnp] dT. ate d (In r) (8) 
(\n p 


> 

where N 5 is the radiance corrected for cloud contamination effects, 
also called the clear column radiance. Since the only constraints 
On ee are that they are known functions of ln p, we can specify, 


(after Ruggles [5] ): 


10 





qT . (V; in p) 
ee d ln p 


~ 


pamilarly, since B[V, T(1n p) He an increasing function of T and 
also of ln p, it satisfied the constraints put on i, in Eq. (5) and 


we can set: 
H, = B[V, TEE p)] 


Substituting these expressions into Eq. (8) produces: 


Ns -( H. F. d Cin p) iC) 


Using this in equation (6) results in a direct expression for ballistic 


tics cerns Of the SLRS radiances. 





v= 3 kx Sarr, N. (10) 
3 Inboe Se 
Now n is the number of SIRS channels used. In practice, the constant 
terms are grouped together and are statistically determined through 
the use of the BIMED 02R stepwise multiple regression routine, (Dixon, 
1 ) available at the W. R. Church Computer Center at the Naval 


Postgraduate School. This produces a final equation of the form: 


Wa-y kxv,[ & b, N, + C (11) 


where b, are the regression coefficients and C, is the regression 


1 
equation constant. The constant terms, g and ln Poo are absorbed by 


the coefficients, b.. 
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B, "D" VALUE APPROACH 

The prackered term nla gies} op (11) represents the regression equation, 
If the ballistic wind, Via is geostrophic, this term can be looked 
upon as a ballistic geopotential height. The regression equation can 
be derived by correlation of SIRS radiances with a dependent sample 
of ballistic geopotential heights. These are obtained by computing 
a weighted average of observed geopotential heights in the vertical. 
The ballistic wind weighting factors are used in the averaging. For 
the dependent sample, deviations from standard, or "D" values of 
height were used for all available mandatory levels from rawinsonde 
soundings. For levels above the maximum height attained by the 
rawinsonde (usually 10 mb) temperature values were extrapolated using 
the standard lapse rate for 45N in January. The hydrostatic 
approximation was then used to calculate "D" values. Linear inter- 
polation was used in the rare instances when an interior mandatory 
level value was missing in a sounding. The following equation 


illustrates the method used for obtaining the dependent sample of 


“Dp values, 


D=S WD (12) 


Here, dD, is the ballistic "D" WenNueg Ue are the ballistic wind weighting 


factors for each mandatory level given by Finke [3], and =o are the - 
actual "D" values of height at each mandatory level. 


This method of calculating V_, then, consists of two steps. 


B? 


First a regression equation is derived by correlation of SIRS radiance 
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data with the dependent sample of "D" values calculated using Eq. 
(12). Then, Eq. (11) is used with the derived regression equation 


to calculate the ballistic wind. 


C. GRADIENT OF RADIANCES APPROACH 


Rearrangement of the order of Eq. (11) yields 
V, = [ee kx 9 N; +C, | (15) 
iz) 


The regression coefficients are represented by C.y and Co is the 


regression equation constant. Here, C. and C, are different than 


2 
b, and C, in Eq. (11). Though equivalent, this is procedurally 
different in evaluation trom Eq. (11). In this case the entire 


expression is a regression equation. The equation is derived by 


; l ANN. 7 = , 
correlation ot values ot ea Witn vi, tne north-south compoonent 
x 


B 
of V,,- Similarly, Ups the east-west component is correlated with 
Van 


The dependent sample of Up and v. is obtained as before 
$ dy B 
from rawinsonde soundings using Kkq. Cie Again, missing interior 
mandatory level values were supplied by linear interpolation. 
Missing upper level values of u were extrapolated by paralleling 
a mean wind profile. Missing v values were extrapolated by simply 
decreasing the already small v component by five per cent with each 


level increment. This was justifiable since v components above 


FO op were Small in all cases. 


D. DIRECT APPROACH 
A purely statistical method suggested by Dr. F. L. Martin was to 


perform direct correlations between the SIRS radiances and the oaliliecie 
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wind components. This produces a very simple regression equation 


for the ballistic winds. 


WV =¢ h. N. +K (14) 


Bese , hy are the regression coefficients and K is the regression 


equation constant. 
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III, SIRS-A EXPERIMENT 


A. DATA AND DATA HANDLING 

The data used in the SIRS-A experiment were obtained from Project 
FAMOS and consisted of radiances, corrected for cloud contamination, 
for the eight SIRS channels for two four-day periods: December 27-30, 
1969, and January 2-5, 1970. The area of data coverage was the 
European continent including the area of the U.S.S.R. west of the 
Urals, and the time period covered was from 0800 GM? to 1200 GM? for 
each day. The dependent ballistic "D" values and winds were computed 
from the 1200 GM? rawinsonde data from the same periods. The 
rawinsonde data did not include any data from France or the British 


Jus aA ellos Jee dba naam 
vu = vw vill ivefors 


Isies. In some caséS5 this was 4 iii 
satellite subtracks lying in this data sparse zone. In addition, data 
was typically quite sparse in the Mediterranean region. 

The 27 degree longitudinal spacing of Nimbus III SIRS data 
precludes the calculation of both u and v components of ballistic 


winds by the geostrophic approximation. Instead, V. the wind 


BN’ 


component normal to the satellite subtrack was calculated using 


Bas. (11) and (13) modified as indicated below 


ae [< p: N: +, | (112) 


eles ; a 
Von aie a d: oui on! (13a) 





Here n is the generally northerly direction parallel to the satellite 
subtrack. The equation for the direct approach remained unchanged. 
The "D" values and values of Up and V3 calculated from all 
available rawinsonde data were plotted and each of the fields analyzed. 
Values were then interpolated to SIRS subtrack positions. Since the 
points with SIRS data were irregularly spaced, the data were processed 
using a cubic spline fitting interpolation routine which yielded 
radiance data at evenly spaced points at a 2 degree latitude interval 
along the satellite subtrack. The effects of this interpolation on 


the methods! accuracy are discussed later. A typical Nimbus III 


subtrack and data is shown in Fig. 1. 


Be. "D" VALUE APPROACH 

The initial regression equations used January 2-4 data as a 
dependent data set. Initial experimentation was confined to the 
January data because the missing French and British data caused a 
bad gap in data coverage relative to the positions of SIRS data on 
December 29. 

At first, two dependent samples of "D" values were correlated 
with radiance data, one using interpolated evenly spaced radiance 
data and another using data from the uninterpolated sub-points. 
Results using the two different data sets were Similar and are 
summarized in Table 1. Since the difference in results was slight 
and the evenly spaced data were much easier to handle, experimentation 


proceeded with the interpolated data. 
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SATELLITE 
SUBTRACK 


© UNINTERPOLATED SIRS-A 
DATA POINTS 


@ INTERPOLATED DATA POINTS 


eee PONEN 
VBN COMPO TS 


Fig. 1. Typical Nimpus III Satellite Subtrack Showing 
Uninterpolated and Interpolated SIRS-A Data. 





TABLE I 


Effect of Interpolation on Results of "D" 
Value Correlations with SIRS Radiances 


Comte lation Suendeardsarror OL 
Cociiicient oCavaiwies Wetameat 
Interpolated 
DEnEE, 2915 51.6 Meters 
Uninterpolated 
Data 979 aoe 


Since the standard errors of "D" value determination were quite 
large compared to those found by Elsberry and Martin [2] fora 
dependent sample of simulated radiance data (22.3 meters), the data 
were carefully examined. It was noted that the great majority of 
points which nad large residuais between rawinsonde derivea "yp" 
values and those obtained from the regression equation were located 
in data-poor regions. It was therefore assumed that the large 
residuals were mainly due to inaccuracies in the dependent "D"' 
values. 

In an attempt to reduce the standard errors of the estimated 
"D" values, those data points which had produced residuals greater 
than 50 meters were removed from the sample. For this filtered 
sample, a new correlation was calculated with somewhat improved 
results, but this reduced the size of the data set from n = 90 to 
n= 49. The correlation was increased to .987, and the standard 
error was reduced to 28.3 meters. 

For the December data, the same two correlation runs were 


performed with the following results: for all data points (n = 80) 
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the multiple correlation coefficient was .806 and the standard 
error was 51.8 meters. For the filtered data (n = 53) the correlation 
was .946 and the standard error was 27.6 meters. The four regression 


equations for ballistic "D" value, D., are listed below: 


B? 


December, unfiltered data: 


dD, = 1.2610, = +2450, + ee + foeenhh - DOs Ne 15) 
+ 05-179, + 9.83O0N_ -—- 1.088N, -— 1822.8 
ii 8 
December, filtered data: 
a -— {. I ° $ - , 
dD, 1 1eN, 1 70eN,, + 14 ofa + 8 JME 58 a eNE (46) 
+ 66.715N¢ + 12.54N,, - 1.834N, - 1986.0 


January, uniiltered data: 


dD. = —9-49 10, - 98. 451N, - 42.484N, + pecan - 1-27 0N, 


: 
+ 40.878N, - 10.980N, + 14-980N, + 1792.3 (17) 
January, filtered data: 
Dy = -5-622N, - 51-070N, + 28.321N, + 20.475N, - 50.7461, i 


+ 19. 161N, + 1.928N. + 12.966N, + 1275.7 


if 


A comparison of the magnitudes of the regression coefficients 
indicated the significance of the different SIRS channels varies 
markedly from the December to the January data sets. In December, 
channels 5 and 6 provide the major modification to the constant 
term while in the January case, channels 2, 3, 4, 5, and 6 all had 
contributions of roughly the same order of magnitude. The coefficient 


for channel 5, however, is almost twice as large as the coetficients 


Le. 





for the other channels. Results reported by Smith, et. al. [6] would 
indicate that the principal contribution in both cases came from the 
layer centered around 250-150 mb. In the January case, however, it 
would appear the principal contribution came from a sligntly lower 
level and that the vertical wind shear was less. 

That the regression equations derived for the filtered and 
unfiltered cases were nearly equivalent was demonstrated by comparing 
"D" values calculated from both periods. The comparisons are 


tabulated in Table II. 


TABLE IT 


Comparison of "D" Values Obtained from Filtered 
and Unfiltered Regressions 


Correlation Explained Variance Soa. errer 
Dec. 993 98.6% 8.7 meters 
Jan. 0995 99.0 20.5 


To check this scheme for calculating usable ballistic winds, 
the "D" values calculated using all four of the derived regression 
equations, Eqs. (15) - (18) were used with Eq. (11b) to generate 
ballistic wind values. These were then compared with the rawinsonde 
derived values as summarized in Table III. 


The mean of the absolute values, \V enh 


of the wind components is 
included in the table since both easterlies and westerlies may be 
included in the sample. In December, the winds were quite variable 
and Wa, and Weyl differ greatly. In January the eae wore highly 


zonal and V.,, andjV,,jare equal. 
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Here, V1 represents ballistic winds calculated using the regression 
equation derived from all data points, and V2, those obtained from 
the equation derived from filtered data. The (I) indicates the 
independent samples, either December 30 or January 5. Although the 
filtered and unfiltered dependent equations were used, both were 
applied to the same, unfiltered data set for the test of the independent 
sample. An unfiltered data set was used as a check since it would 
be impossible to perform the filtering process on an independent 
data set if the procedure were to be used operationally. 

It was found that the significance of the channels entering into 
the dependent regression equations, Eqs. (15 - 18), became very slight 
after more than ae channels were used. The correlation coefficient 
was increased by less than 1.0 per cent with the addition of each 
Succeeding channel. In spite of this, regression equations were 
derived using all SIRS channels since the addition of succeeding 
channels reduced the standard error. 

It can be seen that the filtering of the data had little effect 
on the final outcome. In fact, the results were opposite in the 
two cases. For the independent December data, filtering slightly 
improved the results, while it slightly degraded them in January. 

In either case, the difference was insignificant, at most amounting 
to a 6.1 per cent difference in explained variance or a 2.3 knot 
difference in standard error. The standard errors of the wind 
component are not overly large considering the wind is a weighted 
average from the surface to 30 km. Im any case, they are well 


within the accuracy criteria set down by Finke [3]. 
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For both periods, the wind estimates obtained by using the 
geostrophic approximation with SIRS derived "D'' values are better 
than those derived from "D" values obtained from the rawinsondes. 
Geostrophic values of Vi. were calculated using "D" values interpolated 
to the satellite subpoints and compared with the calculated values 
of Vay at those points. A correlation of .725 was obtained for the 
December data and .766 was the correlation for January. An explanation 
of the superior performance of the SIRS techniques is that the "D" 
value fields obtained from SIRS data are tree of the inaccuracies 
inherent in a field obtained by many different sensors, Unlike 
separate rawinsondes, the SIRS instrument tends to make all errors 


in the same direction, thus not significanly affecting gradients. 


C, GRADIENTS OF N. 


This approach is somewhat more direct than the "D" value approach. 


Here, the derivatives + oN » were compared directly with the 
“ 


Son 


calculated values of V5 and the resulting regression equations were 


used to calculate Ve for the independent sample. The regression 


equation derived from the December data is given below as an example. 


Vix = ~26911X, - 11.456x,, + 30./03x_, - 1.548x 


BN 


p) S 


+64.865x, - 25,090x, = 2.122x, + 2.251x, + 2h.he (19) 


2 { 


a oN: 
Here, x. indicates a Sl ° 
This approach was tried with data interpolated to evenly spaced 
grid points and also with the derivatives calculated at the mid-points 
between the actual, uninterpolated satellite subpoints. In view of 


the results obtained during the first phase of the experiment, no 
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data filtering was attempted. The results of this approach are summarized 
in Table IV. Tests of the December 30 and January 5 independent 
samples for the in situ SIRS data were not made in view of the more 
favorable results for the dependent data when subtrack interpolation 
was employed. 

The results here are somewhat mixed. Only for the December 40 
case does the scheme show as much accuracy as the results for the 
"D" value approach shown in Table III. The differences between the 
results of this method and those of the "D" value approach may be 
explained by the greater number of derivatives which were taken in 
the + ou approach, In the latter case, a derivative was 
calculated for each SIRS channel, while in the "D" value approach 
the only derivative calculated was that of the ballistic "D" value, 
The smaller explained variance with the 7 approach propapdly 
resulted from a summation of the roundoff errors in each derivative 
calculation. In this approach, interpolation to obtain uniform 
Spacing of data produced a greater improvement in results than it did 
in the "D" value approach. This is perhaps explainable by the 
very small values of at compared to the magnitudes of 22s we Gh 
such small values, the loss of definition due to the wide spacing of 
some of the uninterpolated data could seriously degrade results. The 
improvement resulting from interpolation is indicative of the cubic 
spline fitting routine's ability to accurately determine the data 


patterns. 
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D. DIRECT APPROACH 
A third method tried was the derivation of a regression equation 


by direct correlation between V., and N,. This approach gave poorer 


B 
results than the previous two methods used, as Table V indicates. 


(See also Tables III and IV.) 


TABLE V 


Comparison of Winds Obtained From Direct 
Approach with Rawinsonde-—Derived Winds (SIRS-—A) 


Corr Expalined Standard Max 
Sample n Coeff Variance ldueigie ng inaiaene 
27-29 Dec 80 . 788 62.2% 12.8kts 34.3kts 
30 Dec 51 ata 5965 18.0 31.1 
2-4 Jan 88 607 36.8 15e6 38.8 
pees 23 SZ 168 oe 21.1 


Except in the case of the December 2/-29 "D" value approach, 
these results are somewhat inferior to those obtained by either of 
the first two methods described. Two factors could account for the 
superior performance in December. First, the value of W,, is smaller 
in December which would tend to make smaller errors. Secondly, 
the patterns obtained on the u and v analyses were more constant 
with time in the December case and thus the regression coefficients 


obtained from the dependent sample would be more applicable to the 


independent data. 
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Be SUMMARY 

Both the "D" value approach and the gradient of Ne approach 
appear capable of providing usable values of W.. The direct 
approach, however, in most cases gave correlaticns too low to be 
considered of value. It is interesting to note that the correlations 
for the independent data sets were higher than for the dependent 
data with the December data. This result was also noted by Elsberry 
and Martin L2] in their use of the same data for ballistic density 
calculations. 

To further check the temporal stability of the techniques, the 
December equations were applied to the January data set and vice 
versa. The resulting values of ballistic wind were then compared to 
the rawinsonde-derived values. Table VI gives a summary of the 
results for the "D" value approach. All runs were performed using 


using equally spaced data. 


TABLE VI 


Results of Temporal Stability Checks 
of the "D" Value Approach 


Corr Explained otd 
sample Coeff Variance lusaclone 
Dec Data 

Jan HEquations 
V1-Vay WES 1.5% 21.0kts 
ea) Dec 30 677 45.9 2250 
Jan Data 
Dec Equations 

zNF 
V1 Van Se ye. 15.1 
V1-V,,.,(I) Jan 5 -370 13.4 tie 


Za 





It can be seen that the temporal she pile: is rather poor. 
The percentages of explained variances were generally so low as 
to be insignificant. It is interesting to note that the best 
performance was for the December 30 data using equations formed from 
the later data set. A possible implication to be drawn from Table 
VI is that the regression equations would have to be continually 
updated using data from a period immediately prior to the day for 
which the independent wind values are to be calculated. 

In the tables of results for the various methods, it may be 
noticed that the sample sizes vary considerably although the samples 
are taken from the same data set. This is the result of the loss 
of end points from the rows of data in the finite differencing 


processes used to calculate derivatives. 
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IV. SIRS-B EXPERIMENT 


A. DATA AND DATA HANDLING 

The SIRS-B data used in the second phase of experimentation were 
taken from the period January 10-17, 1971. Data were selected from 
the latitude belt from 4ON to 75N. As was the case in previous 
experimentation, corrected clear column radiances were used. 
Concurrent rawinsonde data were also obtained, Though it was 
expected to sacrifice some data density to obtain the side-looking 
capability of the SIRS-B, data coverage was still generally poor. 
Data points ~— widely scattered and many areas of good rawinsonde 
coverage were completely devoid of the SIRS-B data. In areas with 
STRS-B coverage, the normal spacing hetween data points was 450 to 
650 kilometers. Another shortcoming of the SIRS-B data was that 
the window channel was inoperative so that only seven channels 
were available. 

To permit the calculation of gradients, the data were chosen 
from areas where several SIRS data points could be found grouped in 
a "data cluster." These clusters consisted of one master data 
point surrounded by three to six other "slave" points. A schematic 
of a data cluster showing a "master" and five "slaves" is shown 
in Fig. 2. The maximum distance between a master and a slave was 
arbitrarily limited to 1000 kilometers. The surface between each 


master-slave pair was assumed to be linear and a plane was fitted 


between each pair. The average slope for all of the master-slave 
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pairs was then taken as the partial enetive with respect to 
the x or y direction for either the "D" value or the radiance 
value at the master point. 

For the period January 10-17, 45 separate data clusters were 
found in areas with good rawinsonde coverage. These clusters 
comprised a total of 194 data points. The first 30 clusters 
covering January 10-14 were used as the dependent sample in all 
experiments and the 15 clusters from January 15-17 were used as 
the independent sample. The side-looking capability of the SIRS-B 
instrument enabled calculation of both the u and v components of 
the ballistic wind. Eqs. (11) and (13) were therefore broken 
down into their appropriate components. Otherwise, computational 
procedures used were essentially the same as those used in the 


SLRS-A experiment. 


bye "D* VALUE APPROACH 

The regression equation for calculation of ballistic "D" 
values was obtained by correlation of rawinsonde-derived "D" 
values with the radiance data for all 130 data points in the 


dependent sample. The resulting equation is given as Eq. (20). 


D, = -10.82N, + 15.23N, - 47-41 a5 95+ 96N7 - 90. 49Ng 


é 4 


= es (20) 


Note that in this equation, channels 3 and 5 were not statistically 


significant. 
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The correlation coefficient between "D" aes and radiances 
for the dependent sample was .906, indicating 81.2 per cent of the 
variance was explained. The standard error was 119.8 meters. Though 
this standard error was quite large, no data filtering was attempted 
due to the already small sample size and the negligible improvements 
with filtering in the SIRS-A experiment. 

"D't values obtained from the application of Eq. (20) to the data 
points in the indpendent sample were next used in Eq. Gay to obtain 
u and v components of ballistic wind. These SIRS=-B derived ballistic 
winds were then correlated with rawinsonde ballistic winds to check 
the accuracy of this method. The results of these correlations for 
the independent sample are tabulated in Table VII. 

It can be seen that the "D" value approach does quite well at 
determining the v compenent of W ... in fact, the results are 
comparable to those for the one-day January dependent sample as noted 
in Table III. The poorer results with the u component may have 
resulted from the relatively small population of sample clusters. 

The method used to calculate gradients, when used over the wide 
spacings involved in a typical cluster might also contribute to 


the low correlations. 


C. GRADIENTS OF N. 
This technique, using Eq. (13), was carried out in essentially 


the same manner as the experiment with SIRS-A data. In this case, 


ON, 


} 
however, the u and v components were correlated with values Te Sa 


LON; - 
and -= [ut respectively. 
$ ax 
The resulting equations for these components, Up and Vz are given 


below. 
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CS 152.2x, - 2 + pAe 1X) o 197- 6x, (1) 


695 « TX + 2920 1X, + 14.74 


Vp, = -133.4y, - 121.0y; + 144.4y), - 146.0y, 


(22) 
+568.8y, - 129.67, + 5595. 


| 1 ON; 'ON; 
As before, the x, represent fi oy and the 4 represent ; a » The 
results using these equations to calculate U3 and Vp are given in 


Table VIII. Again, the regression technique gave relatively poor 
results for the calculation of Up in the independent sample. The 
high correlations eas Up in the dependent sample indicate that the 
low correlation for the independent sample may be due to data short- 
comings as previously indicated. Comparison with Table IV shows 
that the resuits for Vp, are better tnan for the January SIRS-A 

case. As was noted in the SIRS-A experiment, this method explained 
a somewhat lower percentage of the variance than did the "D" value 


approach, 


Dee LRECT APPROACH 
Direct correlations of the seven SIRS channel radiances with 


both Up and v.. were calculated and the resulting regression 


B 
equations were used as with the SIRS-A data. The results are 
tabulated in Table IX. Note that in this approach the sample sizes 
are much larger than in either of the other methods. 

This method also specifies Up much less accurately than it does 


Vv. The results for both components are not as good as those obtained 


ine 
with the SIRS-A data. (See Table V.) 
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TABLE IX 


Comparison of Ballistic Winds Derived 
by Direct Approach with Rawinsonde—Derived 
Winds (SIRS-B) 


Dependent Corr Explained Std Max 
sample n Coeff Variance Error [Haanatow e 
u component 130 . 288 8. 3% 21kts -47kts 
v components 130 on) 28.6 20 —52 
Independent 
Sample 
U component 64 e 502 Sle 1 } 7 41 
Vv component 64 AG 2074 1S 40 


Because of the wide geographical distribution of the SIRS-B data 
the data sample is much more heterogeneous. The results of this 


test seem decisively insignificant. 


Ke. SUMMARY 

in neither the Up nor v, case were the results of the SIRS-B 
experiments as good as those obtained with the December SIRS-A data. 
For the v component the results were considered quite good considering 
the lack of window channel data and the wide spatial and temporal 
distribution of the data. Since all three approaches gave significantly 
poorer results for the u components, it is thought that a major portion 
Cue Mewatittewty must ae either in the nature of the u field or 
in the small size of the independent sample. 


When correlations between the rawinsonde-derived ballistic winds 


and winds geostrophically calculated from analyzed rawinsonde-derived 
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"D" values were computed, the results were ere For the entire 
test period (n = 45) the correlation coefficient for Up was ~848 and 
for Vz was -815. Making a similar test of correlation using only the 
data of Jan. 15-17 (n = 15), however, the correlation coefficient for 
Up was -Olce and for Vy was .441. These anomalous results indicate 
the variability of the data in the period of the SIRS-B experiment 


and the need for more closely spaced SIRS data points in future 


experimentation. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Both the "D'" value and gradient of radiance approaches were 
found to be capable of specifying the ballistic wind within the 
error limits set down by Finke [3]. Generally, the "D" value 
approach explained a higher percentage of the variance and gave 
lower standard errors. 

Contributing factors to the generally poor results for the u 
components upon independent data testing in the SIRS=-B experiment 
may have been one or a combination of the following. The nature 
of the Up field could be such that it is not well specified by a 
Weegee. ©erress. On eqm@uilon. It was moved that the rogressicn 
equations usually failed to calculate any negative values of Ups itjal 
only one of the two cases where negative values were calculated 
was the Up obtained from rawinsondes also negative. Also, the 
previously exhibited poor temporal stability of the regression 
equations could be more proneunced for Up in this data sample. 

The direct approach is not thought to be capable of claculating 


ballistic winds accurately enough to be useful. 


B. RECOMMENDATIONS 

It is recommended that experimentation be continued with a new 
SIRS-B data set having a more even data distribution around the 
northern hemisphere. Better data coverage would enable the use of 


shorter periods of time to derive the equations which would help 
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counteract the temporal stability problem. Additionally, rocketsonde 
data should be used in addition to rawinsonde data to more accurately 
define the motion and mass fields at upper levels. This, however, 
tends to limit data to the northern hemisphere. 

Experimentation with a long data record applied to homogeneous 
latitude should be conducted to determine the optimum time period 
for determining the regression equations, An optimum record length 
should be found which balances sample size and representativeness of 
the early portions of the record. 

In further SIRS-B experimentation, a more sophisticated method 
for calculation of gradients should be utilized. One method 
suggested by Kung La] is to find the least squares approximation to 


the solution of a system of equations of the following form: 


BQ: = ——— Bx: + 2Q AYi C25) 


Here Q is any scalar variable, and Ax and 4y are the eastward 


th pairing of the master data 


and northward distances between the i 
point and a slave. The unknown partial derivatives may be calculated 
if i 2 2 in each master-slave cluster. The greater the degree of 
overspecification, the better should be the least squares fit. 

Another method for obtaining gradients might be the description of 
the scalar surface as a polynomial by the of a regression equation 
in x and y. Derivatives could readily be calculated from the poly- 


nomial. This method would be time consuming since a separate 


regression equation would have to be generated for each set of points. 
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